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Abstract 

We present a calculation of Wigner distributions for gluons in light-front dressed quark model. We 
calculate the kinetic and canonical gluon orbital angular momentum and spin-orbit correlation of the 
gluons in this model. 
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I. INTRODUCTION 


To have a complete understanding of matter at subatomic level, it is important to understand 
the nucleon spin structure. This means an understanding of how the spin (|) of the nucleon 
is shared by the quarks and gluons in the nucleon and what is the contribution of their orbital 
angular momentum (0AM). Earlier it was believed that all nucleon spin is carried by quarks. 
EMC experiment showed that contribution of quark spin to nucleon spin is very small. So, 
an important question is from where does the missing (remaining) angular momentum comes 
from? As the nucleon is made up of quarks and gluons, it is natural to expect that the missing 
angular momentum come either from gluon spin or from quark or gluon 0AM. This is expressed 
by spin sum rule [1] 


^ = +^' 5 + ^ 

Q Q Gluon OAM 

quark OAM 

Here | Ylq quark and gluon spin angular momentum respectively. The above 

sum rule is called canonical spin sum rule. Except quark intrinsic part, other terms depend 
upon specihc gauge choice. 

Recently, Chen et al [2] introduced a gauge invariant extension (GIE) which is basically a 
prescription to hnd a manifestly gauge invariant quantity that coincides with a gauge non¬ 
invariant quantity in a particular gauge. In this way, one can extend the validity of a physical 
interpretation of the canonical OAM and in light front gauge to any other gauge. There 
is another decomposition called kinetic decomposition of nucleon spin [3] 

^ + L'? + 

where 5'^ and are quark spin and kinetic orbital angular momentum respectively. is 
total gluon contribution to the angular momentum of nucleon. Later, Wakamatsu |1] further 
separated the into a gluon orbital part (L^) and an intrinsic part (S'®) using a prescription 
similar to [2]. Which of the two dehnitions of the OAM is ’’physical” is a matter of intense 
debate. According to the current understanding the difference between kinetic and canonical 
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0AM is the choice of the Wilson line in the dehnition of non-local quark operator: a staple type 
gauge link gives canonical 0AM whereas a straight line gauge gives kinetic 0AM. Interesting 
physical interpretation of both type of 0AM is given in EE] 

Recently it has been shown that the Wigner distribution can provide useful information on 
the spin and angular momentum correlation of quarks and gluons in the nucleon. Wigner dis¬ 
tributions [7j are quasi-probabilistic distribution in which both position and momentum space 
information are encoded. They are directly related to generalized parton correlation functions 
(GPCFs) 0. GPCFs are fully unintegrated off-diagonal quark-quark correlator and contain 
maximum amount of information about the structure of nucleon. If we integrate GPGFs over 
light cone energy {k~), we get generalized transverse momentum dependent parton distribu¬ 
tion functions (GTMDs). These GTMDs are Fourier transform of Wigner distributions and 
vice versa. 

In fact, the GTMDs are directly related to generalized parton distributions (GPDs) |3] and 
transverse momentum dependent parton distribution functions (TMDs) |H], both of which have 
been found to give very useful information about the structure and spin of the nucleon in terms 
of quarks and gluons. We can take two-dimensional Fourier transform of GPDs with respect 
to the momentum transfer in the transverse direction to get impact parameter dependent 
parton distribution functions (IPDPDFs) [10]. These give the correlation in transverse position 
and longitudinal momentum for different quark and target polarization. On the other hand, 
TMDs provide momentum space information, and also the strength of spin-orbit and spin-spin 
correlation. It has been shown that both of these distributions are linked to GTMDs. So we 
can consider the GTMDs, or equivalently Wigner distribution as mother distributions. 

Wigner distributions described above, are joint position and momentum space distributions 
of quarks and gluons in the nucleon. Due to uncertainty principle, they are not positive dehnite 
and do not have probabilistic interpretation. However, integration of Wigner distributions over 
one or more variables relates them to measurable quantities m- Here model calculation of 
Wigner functions are important as these calculations play a signihcant role in revealing what 
kind of information they can provide about the correlations of quarks and gluons inside the 
nucleon. Also such calculations are useful to check various relations among the GTMDs and the 
TMDs and GPDs. In fact some relations have been found to hold only in certain class of models 
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[12]. In this work, we calculate the Wigner distributions in light front Hamiltonian approach 
[131 . This approach gives a intuitive picture of deep inelastic scattering processes, as it is based 
on held theory but keeps close contact with parton ideas [Tl|. Here the partons i.e quarks and 
gluons are non-collinear, massive and they also interact. The target state is expanded in Fock 
space in terms of multi-parton light front wave functions (LFWFs). The advantage of using 
light front (inhnite momentum frame) formalism is that such wave functions are boost invariant, 
so one can work with a hnite number of constituents of the nucleon and this picture is invariant 
under Lorentz boost [13]. In order to obtain the LFWFs of the nucleon, one needs a model 
light front Hamiltonian. However many useful information can be obtained if one replaces the 
bound state by a simple spin-^ composite relativistic state like a quark at one loop dressed 
with a gluon [MllIT]. In our previous work [TH], we studied the Wigner distribution of quarks 
for a simple relativistic spin-^ composite system, namely for a quark dressed with a gluon, 
using light front Hamiltonian perturbation theory. Here we calculate the Wigner distribution 
for gluons in the same model. This calculation is useful as in most phenomenological models 
commonly studied, gluonic degrees of freedom are not present HU, and a study of the gluon 
spin and 0AM is not possible there. 

II. WIGNER DISTRIBUTIONS 

The Wigner distribution for gluons can be dehned as [I9] 




( 1 ) 


where Aj_ is the transverse momentum transfer from the target state and b± is 2 dimensional 
vector in impact parameter space conjugate to Aj^. We calculate Eq. for F*-^ = and 
ri = -ief. 



We have. 


F+* = a+W - a*A+ + gr^^A+A\ 


4 








The gauge field can be written as [H], for i = 1, 2, 


sr^ f dk^cPkx 

( 2 )= 5 :y 

A 

We choose the light front gauge, A'^ = 0, and take the gauge link to be unity. 

In our previous work [18] we calculated the quark Wigner distributions for a quark state dressed 
with a gluon. In this work we investigate the gluon Wigner distribution in the same model 
using light-front Hamiltonian perturbation theory. The state can be expanded in Fock space 
in terms of multi-parton light-front wave functions (LFWFs) as [20] 


,{k)ax{k)e 


-t-"' + eX(k)al(k)e 


^k.Z 




p' ,p±,(J) = <^>'^{p)bl{p)\0) + Mpi] j [dp2] Vl67r3p+<^^(p - Pi - P 2 ) 

(7l(T2 

(p; Pi^P 2 )bi. ipi)ai. (P 2 )|0); (2) 


where [dp] = and ai and a 2 are the helicities of the quark and gluon respectively. 

The LFWFs (d’'^(p) and *^^^ 0 - 2 ) appearing in Eq. (j^ are calculated by solving the light-front 
eigenvalue equation in the Hamiltonian approach. <I)'^(p) is the single particle (quark) LFWF 
gives the wave function re-normalization for the quark and ^aia 2 particle (quark- 

gluon) LFWF. ^'^-^^^{p',Pi,P 2 ) gives the probability amplitude to find a bare quark having 
momentum pi and helicity cxi and a bare gluon with momentum p 2 and helicity a 2 in the 
dressed quark. The two particle LFWF is related to the boost invariant LFWF; gj_) = 

-y/^. Here we have used the Jacobi momenta {xi^qi 2 _) : 


pt = Xip'^, qi± = ki± + Xip_i_] 


( 3 ) 


so that = 0- These two-particle LFWFs can be calculated perturbatively as 


v]/'^ 


,(x,g±) = j 


m 

(1± 


2 _ m‘^+{q±p 

X 


iuP 

1—x 




--T\U 


\/r 


X 


{a±.q±)a± ima±{l - x) 


Xa{e±a2T- 


( 4 ) 


1 — x X X 

We use the two component formalism m X is the two component fermion spinor. T“ are the 
usual color SU{3) matrices, m is the mass of the quark and €±^2 is the polarization vector of 
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the gluon. 


The gluon-gluon correlator in Eq. Q for a quark state dressed with a gluon can be expressed 
in terms of the overlap of two-particle LFWFs. The single particle sector of the Fock space 
expansion contributes only at a; = 1 and we exclude this. 

For T*-^ = <5®-^ we get 


wr'{x,k^,b^) = - 


d^A 


^ J 2(27r)2 






0-10-2 


yX, q±} 


1 *1 I 2 *2 

^C 72 ^Ai ' 


'02'^Ai 


; ( 5 ) 


and for T®-^ = —ie'l we get 


wr'{x,k^,b^) = -t / 


d^A 


<Tl,(725^1 


2(27r) 


J^g-iAx-bx 




1 *2 _ 2 *1 
^02^Ai ^02^Ai 


;( 6 ) 


where x = (1 — x) and q± = —q±. Jacobi relation for the transverse momenta in the symmetric 
frame is given by q'j_ = k± — ^{1 — x) and q± = k^ + ^{l — x). We represent the gluon Wigner 
distribution as where A and A' are polarization of target state and gluon respectively. 


We consider only longitudinally polarized target state and then we have four gluon Wigner 
distributions as follows, in a manner similar to quark Wigner distributions m 
Wigner distribution of unpolarized gluon in unpolarized target state as 

= wr{x, k^, b^) + wt\x, k^, b^); (7) 

Wigner distribution corresponding to the distortion due to longitudinal polarization of the 
target as 

= Wl\x, k^, b±) - Wi\x, k^, b±); ( 8 ) 

Wigner distribution corresponding to the distortion due to longitudinal polarization of the 
gluons as 

= W^{x,k^,b^) + W^{x,k^,b^); (9) 
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and the Wigner distribution corresponding to the correlation due to longitudinal polarization 
of the target state and the gluons 

= W^{x, k±, b±) - fcx, bi_). (10) 

The hnal expressions for these four gluon Wigner distributions are given by, using the 
two-particle LFWFs 


r r nnctA h \ r—^\{Qxq'\){x^ —+ 2) + rn^x^] ■ 

= N dAj - 1 ; (11) 


'D{q^)D{q'^ 


x^{x — 1)^ 


W^^{x,k±,b±) = N I dA^ I dAy 


sin(A_L ■ 6 _l) r 4(2 - x){q 2 q'i - qiq' 2 ) 

-D(<?±)-D(g(|_) I- x‘^{x — 1)2 


( 12 ) 


W^^{x,k^,bA_) = N 


dAx 


dA, 


sin(Aj_ ■ 6j_) |■4(a; 
'D{q^)D{q'^ 


2 x + 2){q2q[ - qiq'2y _ 

x‘^{x — 1)^ -I ’ 


(13) 


W^^{x,k^,bi_) = N dA^ dA, 


cos(A^ ■ b^) r4((g±?i)(2 -x)+ mV) 


D{qx)D{q'^ 


x‘^{x — 1)2 


(14) 


where A^, Ay are x, y component of A± and 


D{ki_) = — 


m? + {k^y (fc±)' 


1 — X 


X 


N = 


2(27r) 


III. GLUON GTMDS AND ORBITAL ANGULAR MOMENTUM 

In order to calculate the gluon GTMDs we use the parametrization as shown in [22] where the 
authors have shown that the correlators like in Eq. ([^ can in general be written as 




2(27r)3 


(p',A' I O(^) I p,A) 


u{p\ A')M^u{p, A); 


(15) 
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FIG. 1: (Color online) 3D plots of the Wigner distributions . Plots (a) and (b) are in b space 
with k± = 0.4 GeV. Plots (c) and (d) are in k space with 6 _l = 0.4 GeV“^. Plots (e) and (f) are in 
mixed space where kx and by are integrated. All the plots on the left panel (a,c,e) are for /S.max = 1-0 
GeV. Plots on the right panel (b,d,f) are for A^aa; = 5.0 GeV. For all the plots we kept m = 0.33 
GeV, integrated out the x variable and we took k_\_ = kj and b± = bj. 
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FIG. 2; (Color online) 3D plots of the Wigner distributions . Plots (a) and (b) are in h space 
with fe_L = 0.4 GeV. Plots (c) and (d) are in k space with 6 _l = 0.4 GeV“^. Plots (e) and (f) are in 
mixed space where kx and hy are integrated. All the plots on the left panel (a,c,e) are for /S.max = 1-0 
GeV. Plots on the right panel (b,d,f) are for A^aa; = 5.0 GeV. For all the plots we kept m = 0.33 
GeV, integrated out the x variable and we took k_\_ = kj and b± = bj. 
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FIG. 3; (Color online) 3D plots of the Wigner distributions Plots (a) and (b) are in h space 

with fe_L = 0.4 GeV. Plots (c) and (d) are in k space with 6 _l = 0.4 GeV“^. Plots (e) and (f) are in 
mixed space where kx and hy are integrated. All the plots on the left panel (a,c,e) are for /S.max = 1-0 
GeV. Plots on the right panel (b,d,f) are for A^aa; = 5.0 GeV. For all the plots we kept m = 0.33 
GeV, integrated out the x variable and we took k_\_ = kj and b± = bj. 
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FIG. 4: (Color online) 3D plots of the Wigner distributions . Plots (a) and (b) are in b space 
with /c_L = 0.4 GeV. Plots (c) and (d) are in k space with 6 _l = 0.4 GeV“^. Plots (e) and (f) are in 
mixed space where kx and by are integrated. All the plots on the left panel (a,c,e) are for /S.max = 1-0 
GeV. Plots on the right panel (b,d,f) are for A^aa; = 5.0 GeV. For all the plots we kept m = 0.33 
GeV, integrated out the x variable and we took k_\_ = kj and b± = bj. 
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where 0(z) stands for the relevant quark or gluon operators and stands for the matrix in 
Dirac space with O determined by the corresponding quark or gluon operator. The amplitude 


shown in Eq. (15) takes the following generic structure when the momentum transfer is purely 
in the transverse direction : 


W. 


ASz,cp 


A', A 


u(p', A) 

TP+ ^ 


(16) 


where cp is the parity coefficient of the partonic operator and ASz is the spin flip number 
given by AS^ = A' — A + AL^ such that A(A') is the initial (final) parton light front helicity 
and ALz is the eigenvalue of the operator ALz = Lz — L'z- Also for twist-2 partonic operators 
we get ALz = 0 and hence ASz = X' — X. 

The gluon operators appearing in Eq. ([^ corresponds to the case when ASz = 0 and cp = ±1 
as shown in eq (3.12), eq (3.42) and eq (3.43) of Eq. ( p)^ . So the relevant parameterization of 
the gluon GTMDs which correspond to T*-^ = 6"^^ in Eq. ([^ is: 


M°’+ 


/ M 




+ 75 


ie 


T 


M2 


cOj+ 

^t,ih 





p0,+ 


+ 


pO,+V . 
M 


where 


ei^ = = -hi. 


(17) 


t-|-l is defined as the twist of the operator in [22], so for twist 2 we take f = 1. 

Comparing and solving Eq. ([^ and Eq. (16), we get the following expression for the gluon 
GTMDs: 
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c0,+ _ 

-2N 

m?x^ + 

— 

2a;+ 2) 

Ai(a;2-2x + 2) 

D{q±)D{q'^) 

x^{x — 1)^ 


Ax^ 

o' 

2N 

1 

to 

1 

_1 




D{q±)D{q'^) 

(1 — x)x‘^ 




p0)+ _ 

■^l,2a 

2N 





D{qx)D{q'^) 

x(fc2Ai - fciA2) 



p0,+ _ 

^l,ib ~ 

2N 




D{qx)D{q'^) 

a;(fc2Ai - fciA2) 




( 18 ) 


with that 

in [H] 

q0,+ 

^l,ia ~ 

p9 . 
-^1,1) 

o' P 

jpg . 

_ 

^ 1,2a 

F 2 ; 

p0:+ _ 

^lAb ~ 

K 

0 


+ Fl3- 

So for the gluon case we get Ff^ as shown below and this relation agrees with that in 

m ^(2 — x) 


(19) 


F9 — 
-^1,4 ~ 


N 

D{q±)D{q'^ 


(1 — x)x‘^ 


( 20 ) 


From this, we can calculate the gluon canonical 0AM since the canonical 0AM is related to 
the GTMD Fi 4 as follows, similar to quarks [HI |21l [2S] 

II = - [ dx(fk±^Ff^. ( 21 ) 
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This gives, 


where, 


ll = -N 


dx{l — x){2 — x) 


Ii — m‘^x‘^12 


( 22 ) 


Ii = 

h = 


cPk I 


m?x‘^ _l_ (fc _|_)2 

(Pk± 

m?x‘^ + {k±)‘- 


= Tilog 


Q'^ + m'^x^ 

g? _|_ 77^2^2 

TT 


{rn?x‘^) ’ 


Here Q and jj, are the upper and lower limits of the k±_ integration respectively. 


For the case when AS'^ = 0 and cp = —1 which corresponds to F*-^ = —ie^ in Eq. ([^ the 
relevant gluon parameterization is given by |22j : 


M°’- = 


/ M \ 



+ 75 


■ krAr 

M 2 





f^T pO- 

\M 


+ 


~W 



(23) 


Again by solving Eq. (23) and Eq. 0 we get the corresponding GTMDs at twist two. 


qO,— 

N 

Ak‘j{x — 2) 

Ai 

D{qp)D{q'^) 

2(a; — l)2a;2 ' [x - 

- l)2a;2 

oO,- _ 
^l,ib ~ 

N 

2(a;^ — 2a; + 2)m^ 


D{aAo{q'^) 

x^(l — x) 


pO- _ 

N 

4(A;i.A i)m^ 

- 

D{q±)D{q'J 

a;(a; - l)(/i; 2 Ai - A;iA 2 ) 

pO- _ 

N 

4A;^m^ 

- 

D{qp)D{q'^) 

x{x - 1)(/i;iA2 - A;2Ai) 


+ A\{x‘^ — 4a; + 5) 
2(a; — l)2a; 


( 24 ) 
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The relation of these GTMDs with that in [S] can be written as : 


cOi~ _ 


1,4’ 


qO,- _ 
pO,- _ 
pO,- _ 


_r^9 . 

^ 1 , 1 ’ 

2m2G?^2 - ^iGh 

2m? 

2m^ 


(25) 


The spin-orbit correlation factor for the gluons Gf can be defined in terms of the GTMD G\ ^ 
as follows [23], similar to quark case [26] 


Cl = / dx(fkx^G\.. 

/ 


(26) 


The GTMD G^ i calculated using Eq. (25) agrees with that 


m 


GU = 


N 


D{q±)D{q'^ 


2(a;^ — 2a; -f 2)m^ 
x3(l — x) 


(27) 


So the spin-orbit correlation factor for the gluons in the dressed quark model is given by : 


C? = -N / dx 


2(x2 -2a; + 2)(l -x) 


X 


Ii — 


(28) 


The kinetic 0AM for the gluons can be calculated using the sum rule for the gluon GPD’s 
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The gluon GPDs in the above relation can be related to the GTMDs as follows:- 



[ d^kxFl,; 

( 29 ) 

E^{x, 0 , t) = 

(fk_i_ 


; ( 30 ) 

H>^{x,0,t) = 

1 d^kxGl^. 

( 31 ) 


Using the above relation and the gluon GTMD calculated above, we can write the kinetic gluon 
0AM in the dressed quark model as 

L’i = ^ J dx{ - + 27r(l - x)|; (32) 

where, 


Q'^ + 

jj2 _|_ ’ 

/(x) = 2x^ — 3x + 2. 

Unlike for the quarks [IH], canonical gluon 0AM and and spin-orbit correlations are different 
in this model. Note that the GTMDs Fi ^4 and Gi,i depend on the gauge link. But up-to 0{as) 
the result does not depend on the choice of the gauge link [23] . 

IV. NUMERICAL RESULTS 

In Figs. [Tj- we have shown the 3D plots for the Wigner distributions of gluon in the impact 
parameter space {bx-by), momentum space (kx-ky) and also in the mixed space {ky-bx)- Normally 
the upper limit of the Fourier transform should be inhnite. But in our numerical calculation, 
we chose an upper limit of | Aj^ | which we called A^ax- Plots on the left and right column 
are for A^ax = 1 GeV and A^ax = 5 GeV respectively. Dependence of gluon Wigner function 
on Amax is similar as the quark Wigner distributions: the peak of the Wigner distribution 
increases in magnitude as Amax increases. First, second and the third row in all Figs. Hl-i 


/i = 


d?k\ 


'im?x‘^ + {kxY 


= ttIoq 


16 











FIG. 5: (Color online) Plots of 0AM (a) Iz and (b) Lf vs niq (GeV) for different values of Q (GeV). 


correspond to the impact parameter, momentum and mixed space plots respectively. The plots 
in mixed space have probabilistic interpretation since we have integrated out the variable in the 
remaining direction i.e. and by and for the impact parameter and momentum space plots 
these remaining variables are held constant. For all the plots, we have taken mass of target 
state to be 0.33 GeV. Also we integrated over x and divided by the normalization constant N. 
In Fig. [T] , we show the 3D plots for the Wigner distribution of unpolarized gluon in an unpo¬ 
larized target state In Figs. 1(a) and 1(b) we see the variation of in the position 

space. The magnitude of is maximum at center = by = Q) and increases with increase 
in Amax, which is expected from the analytic expression of . In Figs. 1(c) and 1(d) we 
have plotted in the momentum space for b± = bj = 0.4. In momentum space too 

peaks at center {kx = ky = Q) and its magnitude increases with increasing Amax- In Figs. 1(e) 
and 1(f) we have shown the variation of in the mixed space. We observed that is 

maximum for ky = 0. As we move away from ky = 0, it hrst decreases and then increases. 
Hence the probability to hnd a gluon in the target state is maximum near ky = 0. 

It is worth mentioning here that similar plots for the quark Wigner distributions in [TH] are 
rotated through an angle | because there we took A_i_ to be positive only whereas here we took 
A± to be both positive and negative. 
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In Fig. we show the 3D plots for the Wigner distribution of unpolarized gluon in longitudi¬ 
nally polarized target state. In Figs. 2(a) and 2(b) we see how varies in position space. We 
observed dipole structure whose magnitude increases with increase in A^ax- In Figs. 2(c) and 
2(d) we have plotted in the momentum space for a hxed h±_ = bj = 0.4. Again we observe 
a dipole structure but the polarity is flipped when compared to the plots in position space. Also, 
the magnitude of the peak increases with increasing A^ax which is expected from the analytic 
expression of . In Figs. 2(e) and 2(f) we have shown the variation of in the mixed 
space. Here we observe quadrupole structure whose magnitude increases with increase in Amax- 

In Fig. 1^ we show the 3D plots for the Wigner distribution showing the distortion due to the 
longitudinal polarization of the gluons. In Figs. 3(a) and 3(b) we show the variation of 
varies in position space for hxed k± = kj = 0.4. We observe that the behavior is similar to the 
case of showing a dipole structure. In Figs. 3(c) and 3(d) we have plotted in the 
momentum space for bj_ = bj = 0.4. In the momentum space we observe a dipole like structure 
again similar to the case of but here the polarity is not hipped unlike that in when 
compared to the plots in the position space. In the mixed space we again observe a quadruple 
structure with increasing magnitude as A^ax increases. 

In Fig. 1^ we show the 3D plots for the Wigner distribution showing the distortion due to the 
correlation between the longitudinal polarization of the target state and the gluons. In Figs. 
4(a) and 4(b) we see how varies in position space. In the b-space the behavior is similar 
to that shown by and the magnitude increases with increasing A^ax value. In Figs. 4(c) 
and 4(d) we have plotted in the momentum space for b± = bj = 0.4. In momentum space, 
shows a behavior similar to that of and its magnitude increases with increasing 

Amax- In the mixed space again the nature is identical to that shown by . 

In Fig[^ we have plotted the 0AM of gluon with respect to mass of the target state for different 
values of Q where Q and fi are the upper and lower limits of transverse momentum integration 
respectively, /i can be taken to be zero if the quark mass is non-zero. In fact we take taken /i to 
be zero. In Figs. 5(a) and 5(b) we show the canonical and the kinetic gluon 0AM respectively 
as a function of the target mass. We observe that the magnitude of both the 0AM decreases 
with increasing mass of target state. 


18 


V. CONCLUSION 


In this work, we presented a calculation of gluon Wigner distributions for a quark state dressed 
with a gluon. This can be thought of as a simple composite spin 1/2 system having a gluonic 
degree of freedom. We showed the various correlations between the gluon spin and the spin 
of the target. We calculated the gluon kinetic and canonical 0AM and also calculated the 
spin-orbit interaction of the gluons. The kinetic and canonical 0AM of the gluons differ in 
magnitude. In most phenomenological models, there is no gluonic degree of freedom and a 
study of gluonic contribution to the spin and 0AM is not possible in such models. Our simple 
held theoretical model calculations may be considered as a hrst step towards understanding the 
gluon spin and 0AM contribution. 
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